Separability Criterion for all bipartite Gaussian States 
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We provide a necessary and sufficient condition for separability of Gaussian states of bipartite 
systems of arbitrarily many modes. The condition provides an operational criterion since it can be 
checked by simple computation. Moreover, it allows us to find a pure product-state decomposition 
of any given separable Gaussian state. Our criterion is independent of the one based on partial 
transposition, and is strictly stronger. 
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Entanglement is the basic ingredient in the philosoph- 
ical implications of Quantum Theory. It also plays a cru- 
cial role in some fundamental issues of this theory, such as 
decoherence or the measurement process. Furthermore, 
it is the basis of most applications in the field of Quan- 
tum Information. However, despite of its importance, 
the entanglement properties of systems are far from be- 
ing understood. In particular, we do not even know how 
to solve the following question [Q: given two systems A 
and B in a state described by a density operator p, are 
those systems entangled? This question constitutes the 
so-called separability problem, and it represents one of 
the most important theoretical challenges of the emerg- 
ing theory of quantum information. 

During the last few years a significant amount of work 
in the field of quantum information has been devoted to 
the separability problem ||. For the moment, the basic 
tool to study this problem is a linear map called partial 
transposition operation. Introduced in this context by 
Peres || , it provides us with a necessary condition for a 
density operator to be separable (equivalently, not entan- 
gled). This condition turns out to be sufficient as well for 
two particular cases: (a) A and B are two qubits or one 
qubit and one qutrit Q ; (b) A and B are two modes (con- 
tinuous variable systems) in a Gaussian state 1|. Thus, 
in these two cases the separability problem can be fully 
solved. However, for higher dimensional systems as well 
as in the case in which A and B consist of several modes 
in a joint Gaussian state, partial transposition alone does 
not provide a general criterion for separability. In both 
cases, examples of states which despite of being entan- 
gled satisfy the partial transposition criterion have been 
provided || fjj. 

In this Letter we solve the problem of separability for 
Gaussian states of an arbitrary number of modes per site. 
Our method does not rely in any sense on the concept of 
partial transposition, and therefore is diametrically dif- 
ferent from the ones that have been introduced so far to 
study the separability problem Q . It is based on a non- 
linear map / : — > "fN+i between matrices jn which 
reveals whether a state p is entangled state or not. Fur- 
thermore, once p is shown to be separable, our method 



allows to find an explicit decomposition of p as a convex 
combination of product states. 

Let us start by fixing the notation and recalling some 
properties of correlation matrices (CMs). A Gaussian 
state of n modes is completely characterized by a matrix 
7 6 M2n,2n (the set of 2n x 2n matrices), called correla- 
tion matrix whose elements are directly measurable 
quantities. A matrix 7 G Mm,in is a CM if it is real, 
symmetric, and 7 — iJ n > 0. Here [H 
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In the following we will consider two systems A and B, 
composed of n and m modes, respectively, in a Gaussian 
state. The corresponding CM will be written as 



7o 



A Q Co 




(2) 



where A 6 M 2nj 2n and B a E M 2 m,2m ar e CM them- 
selves, Co <E Af 2 „ 2 m and J n ,m = J n ® Jm- In order to 
simplify the notation, when it is clear from the context 
we will not write the subscripts to the matrices J and we 
will not specify the dimensions of the matrices involved 
in our derivations. In it was shown that a CM of the 
form (|J) is separable (i.e., it corresponds to a separable 
state) iff there exist two CMs, ja,b, such that 



70 > 1A © 7-B- 



(3) 



This condition, even though it can be very useful to show 
that some particular states are entangled [M |l(| , cannot 
be directly used in practice to determine whether an ar- 
bitrary state is entangled or not, since there is no way of 
determining ja.b in general. If one can determine them, 
however, then one can automatically construct an explicit 
decomposition of the corresponding density operator as 
a convex combination of product states . 

Below we will present a criterion which allows us to de- 
termine whether a given CM, 70, is separable or not, and 
which allows us to determine a product-state decomposi- 
tion if this is the case. To this aim, we define a sequence 
of matrices {7at}^ =0 of the form (S). The matrix 7jv+i 



2 



can be determined through the discrete map defined as 
follows: (i) if 7at is not a CM then jn+i — 0; (ii) if jn 
is a CM then 

A N +i = B N +i = An- Rc(X n ), (4a) 
C N+l = -lm{X N ), (4b) 

where X N = C N (B N - iJ) -1 C£ @. Note that for N > 
1 we have that An — A N = Bn and Cat = — C N are real 
matrices. The importance of this sequence is that, as we 
will show below, 70 is separable iff jn is a valid separable 
CM. In particular, for some finite number of iterations 
7jv will acquire a form in which it is simple to check that 
it is separable. Furthermore, starting from that CM we 
will be able to construct the CMs ja,b of Eq. (||) for 
the original 70. Now we will present several propositions 
from which the above results will follow. Two lemmas 
are presented in an appendix. 

First we show that if 7at is separable, so is Jn+i- More- 
over, the CMs 7a, s associated to jn [cf. Eq. (g)] allow 
us to construct the corresponding CMs for 7jv+i- 

Proposition 1 If for some CMs, Ja,b, we have jn > 
1A © IB then 7at + i > j A © 7A- 

Proof: We use the equivalence (i)-(iii) of Lemma|l|to 
obtain that Bn — C n (An — 7a) _1 CW > 7s > iJ, where 
the last inequality follows from the fact that 75 is a CM. 
Using the equivalence (ii)-(iii) of Lemma ^ we obtain 
7a < A N - C N (B N -«J) _1 Cjy = A N +i +iC N +i, where 
we have also used the map (fjj) . According to Lemma ||, 
this immediately proves the proposition. ■ 

Now, we show that the converse of Prop. [I] is true. 
That is, if 7_/v+i is separable, so is jn- Apart from that, 
the following proposition exhibits how to construct the 
matrices ja.b [cf. Eq. ([5])] related to jn starting from 
the ones corresponding to 7^+1- 

Proposition 2 If for some CM, ja, we have that 
iN+i >JA®7A then Jn >1a®1b, where 

1B =B N - C N {A N - iaV^I, (5) 

is a CM. 

PROOF: We use Lemma |2| and the map (^) to trans- 
form the inequality 7at+i > 7a ©7A into An — CJj(Bn — 
iJ)~ 1 Cn > Ja- According to the equivalence (ii)-(iii) of 
Lemma [j] this implies that 73 > iJ. Since it is clear from 
its definition (|^), 73 is also real and symmetric, it is a 
CM. On the other hand, using the equivalence (i)-(iii) of 
Lemma ^ we immediately obtain that jn > 7a © 7b- ■ 

Using the fact that for N > 1, An = Bn and the 
symmetry of the corresponding matrix jn we have 

Corollary 1 Under the conditions of Prop. II, if N > 1 

we have thatjN > 7a©7Aj where ja = (7a+7b)/2 > iJ 
is a CM. 



The above propositions imply that 70 is separable iff 
7jv is separable for all N > 0. Thus, if we find some jn 
fulfilling (||) then 70 is separable. Thus, we can establish 
now the main result of this work. 

Theorem 1 (Separability criterion) 

(1) If for some N > 1 we have An iJ then 70 is not 

separable. 

(2) If for some N > 1 we have 

LN=AN-\\C N \\o P i>iJ (6) 

then 70 is separable fl^j . 

PROOF: (1) It follows directly from Prop. |l|; (2) We 
will show that 7jv > Ln®Ln, so that according to Prop. 
H 70 is separable. We have 

T m T _L ( I |Cat| |opl C N \ / 7 x 

1^L N ®L N +^ c t {mwi ) , (7) 

so that we just have to prove that the last matrix is posi- 
tive. But using Lemma |TJ this is equivalent to ||Cjv||o p l ^ 
C n Cn, which is always the case. ■ 

This theorem tells us how to proceed in order to de- 
termine if a CM is separable or not. We just have to 
iterate the map (||) until we find that either An is no 
longer a CM or Ln is a CM. In the first case, we have 
that 70 is not separable, whereas in the second one it is 
separable. If we wish to find a decomposition of the cor- 
responding density operator as a convex set of product 
vectors we simply use the construction given in Corollary 
[I] until N — 1 and then the one of Prop. |2|. This will give 
us the CMs, Ja,b, such that 70 > 7a © 7s, from which 
the decomposition can be easily found M . 

In order to check how fast our method converges we 
have taken families of CMs and applied to them our cri- 
terion. We find that typically with less than 5 iterations 
we are able to decide whether a given CM is entangled 
or not. The most demanding states for the criterion are 
those which lie very close to the border of the set of sep- 
arable states (see Corollary || below) . We challenged the 
criterion by applying it to states close to the border of the 
set of separable states and still the convergence was very 
fast (always below 30 steps) . Figure |l| illustrates this be- 
havior. We have taken n = m — 2 modes, an entangled 
CM 7a of the GHZ form |TJ (Fig. |a) and an entangled 
CM 7b with positive partial transposition (Fig. |l|b). 
We produced two families of CMs as 7a,b(e) = j a ,b + el. 
We have determined e a .b such that the CMs become sep- 
arable. In the figure we see that in both cases, as we ap- 
proach exponentially fast e Q , , the number of steps only 
increases linearly. We have also added, instead of 1 other 
positive projectors with all possible ranks and found the 
same behavior. By taking other initial CMs we also find 
the same results. 
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FIG. 1: Number of steps as a function e for CMs of the form 
y„,.h(e) = 7a, t + el where: (a) 7a taken from Eq. (1) in Rof. 
jHf with r = 1/4, and e a = 0.305774915510(1); (b) 7b taken 
from Eq. (9) in Ref. and e b = 0.0978667902228(4). 



Even though we have tested numerically the rapid con- 
vergence of our method, we still have to prove that, ex- 
cept for a zero measure set, it can decide whether a CM 
is entangled or not after a finite number of steps p4[ . 
We start out by considering the set of separable states. 
Since for them 70 > 7_a © 7b with 7,4.5 > iJ, if we just 
consider those with 7^ > iJ , we will leave out a zero 
measure set. In this case we can show that after a finite 
number of steps these separable states will be detected 
by using our procedure. 

Proposition 3 If jq > ja ©7s with ja > iJ + el, then 
there exists some 

JV < JVb = i (||i4o| |tr-2n) + l, (8) 

for which condition (j^) is fulfilled. 

PROOF: Using Prop. |l| we have that for all N, 

An — iJ > el. (9) 

Thus, < Re(Xjv) = An — Ajj+i- Since all the matrices 
in this expression are positive, taking the trace norm we 
have HAvlltr - ||4v+i||tr = ||Re(Xjv)||tr- Adding both 
sides of this equation from N = to iVo, taking into 



N)\\op 



> 



account that ||...|| tr > ||...|| op , and ||Re(X 
||Cjv+i Hop [since Re(XAr) > ±iIm(Xjv)], we have 

Wo-l 

J2 IICW+lllop < \\MW - PjVolltr < Po||t r - 2fl, 
N=0 

(10) 

where the last inequality is a consequence of the fact that 
An > iJ for all N. Thus, among {C'n}n° = i there must 
be at least one for which ||Cjv||op < e. Thus, An — 
1 1 Cn 1 1 op 1 > An — el > where for the last inequality 
we have used Eq. (^), and therefore, for that particular 
value of N, condition (^) must be fulfilled. ■ 
It is worth stressing that from the proof of Prop. |3] it 
follows directly that if 70 is separable, then the sequence 



7n converges to a fixed point 7^ = A^ © B^,, where 
^4-oo = -Sqq > iJ are CMs. On the other hand, for the 
sake of completeness, we will now show that if 70 is not 
separable, then we can always detect it in a finite number 
of steps. We will use the fact that the CMs of insepa- 
rable Gaussian states form an open set, a fact that can 
be directly inferred from condition ([}]). This means that 
if 70 is inseparable, there always exist some eo > such 
that if e < eo then 70 + el is still inseparable and there- 
fore condition (^|) is never fulfilled. However, if 70 was 
separable, then, according to Prop. ||, 70 + el should ful- 
fill that condition before reaching N = Nq. This can be 
summarized as follows. 

Corollary 2 If 7 is inseparable then there exists some 
e > such that starting out from 70 = 7 + el, condition 
^) is not fulfilled for any N < N = {\\A Q \\ tr - 2n)/e. 

Together, Prop. || and Corollary |^ indicate that 
whether 70 is separable or not, and except for a set of 
zero measure, we will be able to detect it in a finite num- 
ber of steps. However, as mentioned above, according to 
our numerical calculations we see that the process always 
converges very fast and in practice one can directly use 
the method sketched after Theorem 1. 

In conclusion, we have obtained a necessary and suffi- 
cient condition for Gaussian states to be separable. The 
condition provides an operational criterion in that it can 
be easily checked by direct computation. It is also worth 
mentioning that our criterion can be used to study the 
separability properties with respect to bipartite splitting 
of multipartite systems in Gaussian states Jl5[] . Our crite- 
rion is based on a non-linear map, and is more powerful 
than partial transposition. This fact indicates that in 
other situations, like the one in which the systems A and 
B are n and m— level systems with nxm > 6, there might 
exist a more powerful criterion than partial transposition 
to determine whether states are separable or not. This 
problem still remains open. However, the results pre- 
sented here represent a significant step in understanding 
the separability problem, which is one of the most chal- 
lenging problems in the field of quantum information. 
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formationsverarbeitung" ) . 



Appendix 

In this Appendix we present the lemmas which are 
needed in order to prove Props. |l| and ^|. 



4 



Let us consider three real matrices < A = A T £ 
M n ,„, < B = B T e M m , m , and C e M„, m , and 



/ \ £ T C B )= M 1 i M,, . ,„.„+,„ . cl:d 



Lemma 1 TTie following statements are equivalent: 

(i) M > 0. 

(ii) kcr(B) C kcr(C) and A - CB^C 7 > 0. 

(iii) ker(A) C ker(C* T ) and B - C T A- l C > 

Proof: We will just prove the first equivalence since 
the other one is analogous. We use that M > iff for 
any two real vectors a £ R n and b £ M m 



a T Aa + b T Bb + a T Cb + b T C T a > 0. 



(12) 



On the other hand, A ~ CB^C 7 > iff for any a £. 
we have 



i T Aa 



i T CB~ 1 C T a > 0. 



(13) 



(i)=>(ii): We assume (|12|). First, ker(i3) C ker(C) since 
otherwise we could always choose a b £ ker(B) so that 
—a T Cb > a T Aa. Second, if we choose b = —B^C 7 a 
then we obtain (|13|). (ii)=^(i): We now assume (13). 
Then, A = CB^C T + P, where P > 0. Defining 
a = B~ 1 C T a, we have that C T a = Ba (since ker(B) C 
ker(C)), and thus the lhs of @ can be expressed as 
a 7 Pa + (a + b) T B(a + b), which is positive. ■ 

In the derivations of Props. |l| and |^ we have not in- 
cluded explicitly the conditions imposed by the present 
lemma on the kernels of B and C. However, one can eas- 
ily verify that all the problems that may arise from these 
kernels are eliminated by using pseudoinverses instead of 
inverses of matrices |]lTf . 

Let us consider two real matrices A — A T £ M n>n and 
C = -C T £ M nn , and 



M = 



A C 
C T A 



= M £ M 2n ,2n- 



(14) 



Lemma 2 M > iff A + iC > 0. 

Proof: This follows from the observation that M is 
real, and that for any pair of real vectors a,b £ R N we 
have (a - ib)\A + iC)(a - ib) = (a © b) T M(a © 6). ■ 
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